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ABSTRACT 

The bright TeV source HESS J1640-465 is positionally coincident with the young SNR G338.3- 
0.0, and the nearby HESS J1641-463 with TeV gamma-ray emission seems to be closely associated 
with it. Based on the nonlinear diffusion shock acceleration (NLDSA) model, we explore the 
emission from these two TeV sources, the particle diffusion is assumed to be different inside 
and outside the absorbing boundary of the particles accelerated in the SNR shock. The results 
indicate that (1) the GeV to TeV emission from the region of the HESS J1640-465 is produced as 
a result of the particle acceleration inside the SNR G338.3-0.0; and (2) the runaway cosmic-ray 
particles outside the SNR are interacting with nearby dense molecular cloud (MG) at the region 
of the HESS J1641-463, corresponding tt^ decay gamma-ray in proton-proton collision contribute 
to the TeV emission from the HESS J1641-463. Also we investigate the possible X-ray emission in 
molecular cloud from synchrotron procedure by secondary produced through escaped protons 
interaction with the MG. 


Subject headings: radiation mechanisms: non - thermal - gamma rays: theory - supernova remnants 


1. Introduction 

It is generally believed that the particles in su¬ 
pernova remnants (SNRs) can be efficiently accel¬ 
erated to relativistic energy through the diffusive 
shock acceleration (DSA) mechanism. Although 
there is no direct evidence to show that the max¬ 
imum energy of accelerated particles in SNRs can 


reach a so called knee energy ^ 10^® eV ('e.g. lBlasi 


20131) . Fermi observations show that the gamma- 


ray spectra of some SNRs (e.g., IC 443 and W44) 
appear the characteristic pion-decay feature, re¬ 
sulting in direct evidence that cosmic-ray proton s 


are accelerated in SNRs ( Ackermann et al.l 1201 111 . 


Therefore, the observations of high energy y-ray 
emissions from SNRs will help us to understand 
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the particle acceleration mechanisms in SNRs. 

Recently, several TeV gamma-ray sources 
(possibly) associated with SNRs in the Galaxy 
have been detected by the Hi gh Energy Stereo¬ 
scopic S ystem (H.E.S.S.) (e.g. lAbramowski et all 


2ni4allt3.l2015^ , which indicate that a SNR with 


TeV emission can accelerate particles up to higher 
energy compared with a SNR without TeV emis - 
sion (e.g., Puppis A ( Abramowski et aP 2015cll l 
if the TeV emission comes from the 7r° decay via 
proton-proton collision. Here we will focus on the 
origins of TeV gamma-rays from two TeV gamma- 
ray sources HESS J1640-465 and HESS J1641-463, 
the latter is very close to the former and is sur¬ 
rounded by a dense molecular cloud (MC). 

HESS J1640-465 is located at the ambient 
region o f a younger (1 1-3.5 kyr) SNR G338.3- 
0.0 (e .g. Shaver fc GossI 1970l : Whiteoak fc Greei] 
1990), it is an extended source with very high- 


energy gamma-ray emission discovered by the 
H.E.S.S during its survey of the Galactic Plane 
in 2004-2006 (Aharonian et al. 2006). Re¬ 
cent study indicate that very high energy (VHE) 
gamma-ray emission from HESS J1640-465 suf¬ 
ficiently overlaps with the north-western part of 
the SNR shell of G338.3-0.0, moreover this VHE 
gamma-ray spectrum naturally connects with 
the GeV spectrum and has a high energy cut¬ 
off ( Abramowski et HI l2014a|) . Gombined with 
the obse rved data of this source b y the XMM- 
Newton, lAbramowski et al.l (|2014al) pointed out 
that the particles can be accelerated up to tens 
of TeV energies in emission region, and the mor¬ 
phology of GeV-TeV emission imply some possi¬ 
ble evidence for proton acceleration in the SNR 
shell of G338.3-0.0. HESS J1641-463 is posi¬ 
tionally coincident with a older (5-17 kyr) SNR 
G338.5-1- 0.1 and appears a hard spe ctrum in VHE 
domain ( Abramowski et all l2014bll . Because of 
its low brightness and proximity to HESS J1640- 
465, the observed feature of HESS J1641-463 
is confusing, and the obvious X-ray candidate 
is not fou nd as an evident a ssociati on so far. 
Therefore, lAbramowski et id] ( 2014lJ l proposed 
three possible scenarios to explain the potential 
origin of VHE hard spectrum of HESS J1641- 
463 based on a simple analytic method. Here 
we will consider the second scenario proposed 
bv lAbramowski et al. ( 2014bll . i.e., the adjacent 
younger SNR G338.3-0.0 would be a major source 


of cosmic rays since SNR G338.5-1-0.1 is older (5-17 
kyr), the escaping protons from the SNR G338.3- 
0.0 interact with the dense MC in the region of 
the HESS J1640-465, subsequently the high en¬ 
ergy gamma-rays are produced through 7r° decay 
in proton-proton collision. 


In this scenario, we need to know not only the 
particle spectra accelerated inside the SNR but 
also that escaping from the SNR. Eor the particle 
acceleration in the SNR, non-linear diffusive shock 

accelerati on (NLDSA) models have been widely _ 

used (e.g. Kang fc Jonesll2006: Zirakashvili fc Aharoni^ 


2010 [ Cajjrioh et aLl20(M Telezhinskv et al. 12012. 


m 


Eerrand et al.ll2(H4 1. which also involved the es¬ 

cape of high-energy particles (e .g. Ellison_&B^d^ 
I 2 OIII : iTelezhinskv et ah 2012bt Kang et al. 2013lb 


Recen tly, based on the NLDSA model of IZirakashvili fc Aharonian 
( 2 OI 2 II and the assumption that the particle dif¬ 
fusion is different inside and outside the absorb¬ 
ing boundar y of the par t icles accelerated in the 
SNR shock, lYang et ^ ( 2015h investigated the 
escape and propagation of high-energy protons 
near young SNRs. In this paper, we will self- 
consistently explain the high energy gamma-ray 
emissions from these t wo TeV sources in the frame 
of the model given bv lYang et al. ( 2015ll . In Sec¬ 
tion 2, we review briefly the NLDSA model used 
here. We show the model application and the cor¬ 
responding results in Section 3. Finally we give 
our conclusion and discussion in Section 4. 


2. The Model 


2.1. Particle Spectrum Accelerated in a 
SNR 


Zirakashvili fc AharonianI (|2012ll described a 
numerical model of the NLDSA in a SNR. In this 
model, hydrodynamic equations for gas denisity 
p{r,t), gas velocity u{r,t), and gas pressure 
Pg{r,t) (see Eqs. (1) -(3) of their paper), com¬ 
bined with the equation for the quasi-isotropic 
particle momentum distribution N{r,t^p) pre¬ 
senting the transport and acceleration of rela¬ 
tivistic particles by the forward shock and re¬ 
verse shocks (see Eq. (4) of their paper), are 
solved numerically in the spherically symmet¬ 
ric case. Since the pressure gradient dPc/dr 
{Pc = 47r f p'^dpvpN{r,t,p)/3) of accelerated par¬ 
ticles with a momentum p and a velocity v is 
included in the hydrodynamic equations (see Eqs. 
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(2) -(3) of [Zirakashvili fc Aharoniaij (|2012f) '). the 
accelerated particles will react onto the shock 
structure, resulting in a self-regulation of ac¬ 


celera tion efficiency. IZirakashvili fc Aharonian 


(1201211 described the numerical procedure of solv¬ 
ing above equations in detail. 

We now present the relevant physical parame¬ 
ters in this model. Firstly, thermal protons with 
the momenta p = Pf and p = Pb are assumed 
to be injected at the fronts of the forward and 
reverse shocks at r = Rf{t) and r = Rb{t), re¬ 
spectively, where the dimensionless parameters rjf 
and rjb are corresponding injection efficiencies and 
subscript / and b for quantities are represented 
the forward and reverse (backward) shock respec¬ 
tively. Secondly, the accelerated protons propa¬ 
gate inside the SNR through advection and dif¬ 
fusion. On the one hand, the advective veloc¬ 
ity of the protons is w{r,t) = u{r,t) + ^aVa/V^ 
in the isotropic random magnetic field R, where 
Va — B/^/Anp is Alfven velocity and ^a is a pa¬ 
rameter which represents the relative direction of 
fluid velocity to Alfven velocity with = — 1 
in upstream of forward shock, = 1 in up¬ 
stream of reverse shock and ^a = 0 in the other 
region (e.g. the downstream for forward and re¬ 
verse shocks ). On the other hand, the diffusion 
coefficient of protons inside the SNR can be ex¬ 
pressed as Df b = where rjE = 2, 

Db = vpc/iqB is Rohm diffusion coefficient, Pf 
and Pb are the momenta separating proton dif¬ 
fusion in forward and reverse shock regions (t he 
detailed description see lZirakashvili et al.l ( 2014ll l. 
Below Pf = 0.61 TeV/c and Pb = 6.1 TeV/c 
are assumed to reproduce the gamma ray obser¬ 
vations. 

Although the modeling of amplification and 
transport of magnetic held cannot be taken into 
account for consistency, the magnetic held has 
been introduced to play dynamical role in down¬ 
stream regions both reverse and forward shock, the 
pressure and energy hux of magnetic held cannot 
be ignored due to strong magnetic held, because 
they have effects on the dynamical evolution of the 
SNRs shock, therefore the magnetic pressure and 
magnetic energy hux sho uld be introduced to the 
huid dynamical equation (j Zirakashvili et al.ll2ni4[l 
in downstream areas, and the magnetic held spa¬ 
tial distribution depending on radius at r > Rf 


is 



Here po is the gas density of the circumstellar 
medium. The parameter Maj determines the 
value of the amplihed magnetic held strength at 
forward shock, in the shock transition region the 
magnetic held strength is derived via Hugoniot 
condition, i.e. in downstream region the magnetic 
held strength is increased by a factor of , 

where a is the shock compression ratio. In the 
upstream region of the reverse shock we also used 
the equation similar to Eq. o, only replaced Vf 
and Maj with 14 and Maj- In following applica¬ 
tion, we use the adiabatic index of isotropic ran¬ 
dom magnetic held 7m = 4/3, the strength factors 
of magnetic held Maj = 4 and Maj = 8. 


2.2. Particle Spectrum Escaped from a 
SNR 

The NLDSA with the reaction of accelerated 
particles onto the shock structure has been pro¬ 
posed in detail to research the radiation process 
in SNRs ((zirakashvili fc Aharonianll2012 ). Subse¬ 
quently, the processes of accelerating particles in 
shock front and the escape of high energy parti¬ 
cles at a abs orbing boundary have been taken into 
accou nt (e.g. Telezhinskv et al.ll2012bl : Yang et al 
201511 . 


To self-consistently explain VHE emissions 
from HESS J1640-465 and HESS J1641-463, we 
deal with the particle diffusion by dividing the 
simulating domains int o two diffusion reg ions (also 
see the hrst model of Yang et al. I l2015f) . one is 
Bohm-like type in the vicinity of the forward shock 
within a radius of 2Rf (inside SNRs), the other 
is generally considered as Galactic diffusion when 
accelerated particles is far away the Bohm limited 
region (e.g. r > 2Rf is outside SNRs). The prop¬ 
agation of high energy particles outside a SNR is 
considered as Galactic diffusion, and the diffusion 
coefficient can be expressed as 


DgiEp) = X 


E-n 


10 GeV 


cm^s ^ 


( 2 ) 


where Ep is the accelerated proton energy, y is 

the correction factor in the case of slow dif¬ 
fusion around the SNRs, its value (y < 1) is 
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lower than that of av erage Galatical diffu - 
sion case (x ~ 1) (see Gabici et al.l ( 20071 
also s ee t heir argumen t s of iLi & Chen 


(2012) and lOhira et al. (2011)), and 6 
0.3—0.7 is the energy-depe ndent index of the 


diffusion coefficient (e. g. (Berezinskii et al. 


1990 ; Fuiita et al. 2009h . Below we use the 
value of ^ = 0.5 and y = 0.1. 

2.3. Intrinsic 7 -ray emission inside a SNR 

With the NLDSA model, the hydrodynamic 
process and particle propagation can be calcu¬ 
lated by means of combining the evolution of 
magnetic field and particle diffusion. Thus, for 
a given SNR, the spectra of accelerated parti¬ 
cles are produced naturally at any time t . In 
general, high energy gamma-rays are thought to 
originate from two scenarios, one is hadronic sce¬ 
nario, the radiation are produced through the 
proton-proton interaction process, and the other 
is inverse-Compton (IC) process by energetic elec¬ 
trons scattering up soft photons, this also called 
leptonic scenario. In this paper, the seed pho¬ 
tons of IC come from 2.7 K microwave back¬ 
ground light, and we will use the emissiv i ty for - 
mula developed by iBlumenthal fc Gould ( 197Clll . 
the proton-proton interaction process which can 
be calculated by the an alytic approximation given 
bv iKelner et al. ( 200bh . 


2.4. Nonthermal photons in a MC nearby 
a SNR 

The escaped protons outside a SNR during 
their propagations in Galaxy are assumed to col¬ 
lide a nearby MC, so the proton-proton interaction 
becomes an important process producing nonther¬ 
mal photons in the MC. Such a process produce 
neutral pions (tt^) and charged pions (tt^). On 
the one hand, high energy g amma-rays can be 
created through 7r° decay (e.g. Kelner et al. 2006t 
Zhang fc Fang Il2007t iTang et al. 2011 ). On the 
other hand, decay will produce secondary elec¬ 
trons and positrons, and these pairs will emit 
nonthermal photons through synchrotron radia¬ 
tion and inverse Compton scattering. Generally 
a molecular cloud has a density nuc > 10^ cm“^ 
and magnetic field Bmc > 10“®G, if the pro¬ 
duced pairs diffuse in the MC with 
a radius RmGi then the typical diffusive 
time-scale is tdif 3 x 


103(Rmc/ 10 pc)2x-i(E/10 GeV)-'5(BMc/100 


yr. At the same time, the energy loss time 
for electron synchrotron emission is given 
by Avn = 1.3 X 103(E/TeV )-i(RMr/100 /rG)-^ 
yr (see Gabici et al.l^09f) . With typical pa¬ 
rameters Rmg ^ 10 pc, X = 0.1, and <5 = 0.5, 
the ratio of tdif to tgyn is approximately 


tdif/Wn ~ 2(£;/TeV)0-5(RMc/100 Obvi¬ 

ously, for strong magnetic field (> 100 nG), 
the synchrotron loss time for electrons 
with their energies ^ TeV is shorter than 
the diffusion time, the secondary electrons 
would radiate all their energy through syn¬ 
chrotron emission, thus the diffusive trans¬ 
port of electrons within the MC can be 
neglected. 


The density of the MC should satisf y 
some distributions ( Protheroe et al. Il 2008 h . 
the density profiles could affect the distri¬ 
bution of magnetic field, particle diffusion 
and the following non-thermal emissions. 
For simplicity, in this paper, we assume a 
fiat density profile of the MC with an aver¬ 
age value of riMC = 10 ^ cm“^. 

After obtaining the production rate, q{E)±, of 
secondary e^, we can readily obtain the number 
densities of electro ns and positrons per unit energy 
( ■lones et al.l 200^ . N^{E), within the molecular 
cloud: 

J^q±iE')dE' 


N±iE) = 


dE/dt 


(3) 


where dE/dt is the total energy loss rate of elec¬ 
trons at given energy E dne to synchrotron emis¬ 
sion. 

dE/dt = , (4) 

where Ub = is the density of mag¬ 

netic field within the molecular cloud. Then 
the synchrotron emissivity ji, is calculated by 
using formulae in synchro tron radiation theory 
( Rvbicki fc LightmanI Il979l ) . At the same time, 
the secondary also generated IG gamma-ray 
emission, the process and seed photons is the same 
as the corresponding description of Sect. 2.3. 


3. Application 

In this section, we apply the model to explore 
the high energy emission origin of the two TeV 
sources: HESS J1640-465 and HESS J1641-463. 
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It is thought that the particles are accelerated 
inside SNR 338.3-0.0, then high energy gamma- 
rays yielded as a result of 7r° decay in the proton- 
proton interaction process in SNR, these gamma- 
rays could just present at the region of HESS 
J1640-465. In addition, some high energy particles 
can escape out the SNR G338.3-0.0 at the region 
of r > 2Rf, the runaway particles will collide with 
the mater in dense MC by pro ton-proton interac¬ 


tion (al so see the discussion of lAbramowski et al 


( 2014bll . and then some high energy gamma-ray 
could be produced at the area of HESS J1641-463. 


3.1. HESS J1640-465 


HESS J1640-465 is a VHE-ray source observed 
by H.E.S.S., which is positionally coincident with 
the SNR G338.3-0.0 and e vidently is a shell-type 
SNR with 8' diameter (e. g. lAharonian et al. 20061 : 
Whiteoak &: Green 1996 ). Based on 21 cm Hi 
absorption spectra to G338.3-0.0 and adjacent 
Hii r egions, the distance to G338.3-0.0 is 12-13.5 


kpc (jLemiere et al 


20091 ). While the distance is 
limited to a small range well, there is no direct 
measurement of the ambient density, but similar 
to RX J1713.7-3946, the HESS J1640-465 may 
eyolve in a wind bubble of progenitor star at 
present epoch, howeyer some eyidences for ther¬ 
mal radio emission in this SNR s hell also indicate 
the p resence of dense material (jGastelletti et al 
20111 ). The age of SNR G338.3 -0.0 is still ambigu - 


ous, and would be (5-8)ky r ( Slane et al. 2010l) . 


recently lAbramowski et ahl ( 2014a ) re-estimate 
the age would be (1-2) kyr which is significantly 
younger than the estimate of (5-8) kyr. Here we 
adopt an age of 2000 yr, a distance of 13 kpc, and 
an ayerage number density of no = 0.2 cm“^. 


In our calculations, the supernoya explosion en¬ 
ergy Esn = 10®^ erg and the initial ejecta mass 
Mej = 1.4M0 are assumed, the SNR shock prop¬ 
agates in the ambient medium with an ayerage 
hydrogen number density no = 0.2 cm“^, a mag¬ 
netic field strength Bq = 1.0 p,G, a temperature 
T = 10^ K, and the index of ejecta velocity dis¬ 
tribution k = 7. Below we use the injection pa¬ 
rameter rjf = r]t, = 5 X 10“^ for protons, and the 
electron injection with rjf = rji, = 10“^. The value 
of other parameters have been described in above 
section. The following results is obtained at the 
epoch t = 2000 yr. 


We present the radial dependencies of the gas 




Fig. 1.— Top panel: radial dependencies of the 
gas pressure (solid line), the GR pressure (dashed 
line), and the magnetic pressure (dotted line) at 
the epoch t = 2000 yr, RS, FS and GD are the 
positions of the forward shock, reverse shock and 
contact discontinue respectively. Bottom panel: 
spectra of accelerated protons (solid line) and elec¬ 
trons (dashed line)in the position of forward shock 
at the epoch t = 2000 yr. 
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pressure, the accelerated particle pressure, and the 
magnetic pressure in the top panel of Fig. [T1 where 
the accelerated particle, gas, and magnetic pres¬ 
sures at forward shock are 30%, 43%, and 8% of 
the ram pressure. In the same figure we also show 
the forward shock (FS), reverse shock (RS) and 
contact discontinue (CD) positions, where the for¬ 
ward shock position Rf = 15.5 pc is at the upper 
boundary of 4' (about 15 pc) of radius for its ob¬ 
served value. The spectra of accelerated protons 
and electrons at the position of forward shock are 
shown in the bottom panel of Fig. [1] the pro¬ 
tons will be accelerated to run up to maximum 
energy ~ 50 TeV, which appear to be a source 
as significant contribution of Galactic cosmic ray 
flux around the knee. Because of the synchrotron 
losses of accelerated electrons in strong magnetic 
field, the maximum energy of electrons were only 
accelerated up to several TeV. 

With above results, we calculate the spectra 
of high energy y-rays which are produced by 
hadronic 7r° decay and leptonic IC scattering re¬ 
spectively. Our results are shown in Fig. [2] The 
solid line represent the emission produced through 
the accelerated protons from G338.3-1-0.0 inter¬ 
acting with the ambient gas, and IC gamma-rays 
(dashed line) have been generated by accelerated 
electrons scattering off soft photons from 2.7 K 
microwave background radiation. For compari¬ 
son, the observed spectrum of HESS J1640-465 is 
plotted in this figure. It can be seen that the ob¬ 
served GeV-TeV spectrum can be reproduced well 
in this model, and has a hadronic origin. Note 
that the leptonic contribution is not clear due to 
poorly X-ray observation. 


3.2. HESS J1641-463 


The protons which are accelerated up to high 
energy in SNR G338.3-1-0.0 can escape from the 
accelerated site into Galactic diffusion region. We 
calculate the escaped proton spectra for different 
diffusion distances at a given time t — 2000 yr. 
Our results are shown in FigEJ it can be seen that 
the energy flux of accelerated particles is sensitive 
to the diffusion distance, particular in the lower 
energy part of the spectrum. HESS J1641-463 is 
located only 0°.25 ( r ^ 50 pc for 13 kpc of a 
distance to observers) away from HESS J1640-465 


( Lemoine-Gonmard. et al. 201^, below we used 


the particle spectra at the diffusion distance r = 


50 pc to calculate the radiations from HESS J1641- 
463. 

If these escaped protons go into a MG nearby 
the SNR, then they will interact with the the MC’s 
medium. For a typical MG, its number density and 
mass of are umc = 100 cm“^ and Mmc = IO^Mq 
respectively, we calculate the spectra of nonther- 
mal photons which are produced by 7r° decay and 
inverse Gompton scattering of the secondary elec¬ 
trons and positrons and show the results with 
the observed spectrum of HESS J1641-463 in the 
left panel of FiglU In this case, the gamma- 
ray emission mainly comes from decay in the 
proton-proton interaction (thick solid line), and 
inverse Gompton scattering of the secondary e^s 
also has a contribution which may dominate over 
lower energy region of gamma-rays and depends 
on the magnetic held strength (thin lines). How¬ 
ever, our predicted flux of high energy gamma 
rays in Fermi/LAT energy range is lower than 
10“^^ erg cm“^ s“^ in the ^ 100 GeV band, a 
factor 10^ lower than the hux of IFHL J1640.5- 
4634 at those energy band, even we plus the con¬ 
tribution from the secondary leptonic IC emis¬ 
sion. A detection of GeV excess on the position of 
HESS J1641-463 would indicate either a contribu¬ 
tion from an unrelated source or from a different 
component of radiation of the same source. 

To explore the possible X-ray radiation from 
the source of HESS J1641-463, it is assumed that 
the X-rays come from synchrotron emission of sec¬ 
ondary electron/positron pairs. The results are 
shown in the right of FigE) Obviously, the syn¬ 
chrotron emission flux level seems featureless, the 
only result is that the peak frequency is moving 
to high frequency with enlarging magnetic field 
strength. Particularly, if the magnetic field larger 
than 1000 ^Gs in the MG, then the peak energy 
will locate at a few keV region, and the flux of 
X-rays may be detected by NuSTAR. 

4. Summary and Discussion 

In this paper, we have investigated the TeV 
emissions from two TeV gamma-ray sources 
(HESS J1640-465 and HESS J1641-463) and given 
a self-consistent explanation in the modified non¬ 
linear diffusion shock acceleration model. On the 
one hand, the TeV emission from HESS J1640- 
465 is naturally produced through 7r° decay in 
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Fig. 4.— SEDs of HESS J1641-463. Left panel: The high energy emission produced through the escaped 
energetic protons interacting with MC (thick solid line), and the secondary electron IC scattering in different 
magnetic field strength (their values are as same as those shown in right panel), observed data of HESS J1641- 
463 derived from Abramowski et al. (2014b). Right panel: possible secondary electron synchrotron emission 
for different magnetic field strength in MC, the thick solid line represents the sensibility of NuSTAR. 



Fig. 2.— The high energy emission spectra of 
HESS J1640-465. Observed data of HESS J1640- 
465 derived from Abramowski et al. (2014a). The 
solid line and dashed line represent the very high 
energy emission from hadronic tt^ decay and lep- 
tonic 1C scattering respectively. 



Fig. 3.— Escaping particle spectra with S = 0.5 
from SNR are illustrated at the distances of r = 
50 (solid line), 100 (dashed line), and 200 (dotted 
line) pc for a given time t = 2000 yr 
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the proton - proton interaction inside the SNR 
G338.3+0.0 (r < 2i?/) and the inverse Compton 
scattering of the accelerated electrons inside the 
SNR has a negligible contribution (see Fig. [2]). On 
the other hand, high energy protons escaped from 
the SNR G338.3+0.0 undergo the Galactic diffu¬ 
sion and collide with the MG which is positionally 
coincident with HESS J1641-463, so the TeV emis¬ 
sion of the HESS J1641-463 could come from the 
runaway protons interacting the MG at the re¬ 
gion of r > 2Rf, the possible X-ray emission of 
the HESS J1641-463 is predicted as being relevant 
to secondary products of the runaway protons, it 
is worth looking forward that future detailed ob¬ 
servations will provide important information for 
transparent explanation. 


Finally, we would like to point out that we 
have given a self-consistent explanation of the TeV 
emissions from two TeV gamma-ray sources in the 
sce nario which is one of the scenarios proposed 
bv lAbramowski et al. ( 2014bll . but other possible 
scenarios are not excluded. 
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